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ABSTRACT 


A method is described to calculate the angle through which a tele- 
scope mounted on a spinning spacecraft "sees" the sunlit portion of a 
specified central body. The central body is assumed to be a sphere and 
the spacecraft is assumed to make a full revolution at the given time. 
The equations of the terminator plane (the plane separating the sunlit and 
dark portions of the central body) and the tangent plane (the plane se- 
parating the portions of the central body that are visible and non-visible 
from the spacecraft) are derived. 
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SIMULATION OF AN OPTICAL ASPECT SYSTEM 


I. INTRODUCTION 

The Optical Aspect System for the AIMP vehicle has three telescopes 
mounted at 45, 90, and 135 degrees to the spin axis. These telescopes are 
connected to photosensors which are activated when the sunlit portion of the 
earth or moon is "seen" by the telescope. This paper will describe a method 
for calculating the angles through which each of the telescopes "sees" the sun- 
lit portion of the central body. This information, together with knowledge of the 
spin axis-sun angle and spacecraft position is sufficient to determine the inertial 
spin axis orientation. This method requires the vehicle's position,^ , the sun's 
position, p s , the unit spin axis vector, S° , and the radius of the central body. 

The method described assumes an instantaneous rotation of 360 degrees about 
the spin axis at the given time. The central body is considered to be a sphere. 
The coordinate system uses the center of the central body as the origin. 


n. TANGENT AND TERMINATOR PLANES 
A. Tangent Plane Equation 

The tangent plane is the plane of the intersection of the cone from the vehicle, 
tangent to the central body, and the surface of the central body (see Figure 1). 
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Let x be a vector from the center of the central body to a point on its sur- 
face which lies in the tangent^plane. Then A = x - H isa vector in the tangent 
plane (see Figure 1) and P * A = 0. 


I H I = lxl sin 6 - R sin 9 = - 


where R = |x| , the radius of the central body. 


\P 


p • A = p * x - 


R 2 £_ 
\p\ \p\ 


= 0. 


Let ~p = (a, b, c) , and x = (x, y, z). Then 


p • x - p • |R 2 


p • x - R-* = 0 


P 


or ax + by + cz r R = 0, which is the equation of the tangent plane. 


B. Terminator Plane Equation 


The terminator plane is the plane which divides the sunlit and dark portions 
of the central body. It is assumed that this plane passes through the center of 
the central body. 

Let p s - (u, v, w) be the vector from the center of the central body to the 
sun, and x = (x, y, z) be a vector from the center to a point in the terminator 
plane. Then x • p s = 0 and the equation of the terminator plane is ux + vy + 
wz = 0. 


m. TELESCOPE PERPENDICULAR TO SPIN AXIS 
A. Equation of Telescope Plane 

The telescope plane is the plane through the vehicle, perpendicular to the 
spin axis vector S°. Let x = (x, y, z) be a vector to a point in the telescope 
plane, and let S® = (£ , m, nt). Then x - p lies in the telescope plane and 
S ° • (x- p) =0. Letting p = (a, b, c), the equation of the telescope plane be- 
comes £x + my + nz - (£a + mb + nc ) = 0. 
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B. Calculations of Intersections 

In the 90 degree case, the points of intersection of the tangent plane, tele- 
scope plane, and the central body are found. If two such intersections are not 
found, the viewing angle is 0. If the intersections are found, they are tested to 
see if they are sunlit and then the intersections (if any) with the terminator plane 
are calculated. (For mathematical solution, see Appendix). 

The tangent plane solutions are tested to see if they are sunlit. The dark 
solutions are discarded. A solution will be sunlit if, any only if, the vector, R s 
from the center of the central body to the solution and the sun vector ~p s form 
an acute angle. That is R s * p s > 0. 

Terminator plane solutions must be tested for visibility from the vehicle. 
The solution is visible if its distance from the vehicle is less than the length 
of the tangent from the vehicle to the central body. 


C. Viewing Angle 

If R Sl and R S2 are two sunlit, visible solution vectors, the viewing angle 
8 is the angle between the two vectors from the vehicle to the solutions, 


8 = cos -1 


(P - Rs x ) ' (ft - Rs 2 ) 

I p - Rsj I Ip- Rs 2 I 


IV. TELESCOPE NOT PERPENDICULAR TO SPIN AXIS 

If the telescope-spin axis angle is not 90 degrees, the calculations are more 
complicated. First the equation of the telescope cone is calculated, and a test is 
made to see if the telescope cone lies outside the tangent cone. Then the inter- 
sections of the tangent plane, telescope cone, and central body are computed 
from a quartic equation. Next a test is made to eliminate spurious solutions 
introduced by squaring. The good solutions are then tested to see if they are 
sunlit. If there are two solutions, the tangent plane viewing angle is calculated. 

If there are fewer than two solutions, the solutions are saved, and the above 
calculations are repeated with the terminator plane replacing the tangent plane. 

If there are two terminator solutions, the terminator plane angle is calculated. 

If there is one terminator and one tangent solution, the viewing angle is calcu- 
lated from these. If there is one terminator and no tangent solutions or no 
terminator solutions, the cone is tested to see if it is sunlit and the viewing angle 
is set to 0 or 360 degrees accordingly. 
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A. Equation of Telescope Cone 


The telescope cone is the locus of points x = (x, y, z) such that x - p and 
S° make an angle equal to the telescope-spin axis angle <*. Thus 


(x-/3)-S°=|x-/?| | S ° | cos a. 

Again let p = (a, b, c), S° = ( £ , m, n). Then 


£ (x - a) + m (y - b) + n (z - c) = /(x - a) 2 + (y - b) 2 + ( z - c ) 2 + m2 + n2 cos a - 

Squaring both sides, and letting the equation of the cone be 

a t x 2 + a 2 y 2 + a 3 z 2 + a 4 xy +.a 5 xz + a 6 yz + a ? x+a 8 y + a g z + a 1Q = 0, 

we obtain by collecting terms 

= £ 2 - (£ 2 + m 2 + n 2 *) cos 2 a, a 2 = m 2 - (£ 2 + ra 2 + n 2 ) cos 2 a, 
a 3 = n 2 - (£ 2 + m 2 + n 2 ) cos 2 a, a 4 =2£m, a 5 =2£n, a 6 = 2mn, 

a ? = - 2 £ (a£ + bm + cn) + 2a (£ 2 + m 2 + n 2 ) cos 2 a 
a g = - 2m (a£ + bm + cn) + 2b (£ 2 + m 2 + n 2 ) cos 2 a 
a 9 = - 2n(a£ + bm + cn) + 2c (£ 2 +m 2 +n 2 ) cos 2 a 

a lo = (a£ + bm + cn) 2 - (a 2 + b 2 + c 2 ) (£ 2 + m 2 + n 2 ) cos 2 a. 

B. Calculations of Intersections 

The points desired are the points of intersection of the telescope cone, 
tangent or terminator plane as the case may be, and the central body. The 
equations of these surfaces have been obtained above. For the mathematical 
solution, see the appendix. 
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C . Tests on Solutions 


For a correct solution the vector from the vehicle to the solution must make 
the specified angle, a , with the spin axis, 


(R s -p) -s° 
!* s -p\ 


- cos a 


The sunlit and terminator plane solution visibility tests are the same as in 
the perpendicular case. 


D. Viewing Angle (Figure 2) 
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If there are two sunlit, visible solutions, take the vectors X x , and x 2 from 
the vehicle to the solutions and project them parallel to the spin axis into a plane 
perpendicular to the spin axis. The angle between the two projections, Px t and 
Px 2 , is the true rotation angle of the vehicle and hence the viewing angle (or, 360 - 
the angle). If fi = v/2 - a, and considering subscripts fixed, |y| = |x| sin (3 
and Px = x - |y \ S° or Px = x - |x | sin fi S° . The viewing angle is either 
cos* 1 (Pxj-Pxjjor 360 - cos -1 (Px x * PX 2 ). To decide which, calculate the mid- 
point of the line joining the endpoints of Px x and Px^, and project it parallel to 
the spin axis onto the telescope cone. The angle isPx x • Px 2 if the vector from 
the vehicle through the point just found is sunlit (see Figure 3) . That is 


Px t +Px 2 


|Q 3 I = iQj tan/3 

Q 2 ± |Q 3 | S° 


Q 2 is tested the same as in testing for a sunlit cone in (E) below. 



Figure 3— Midpoint Test 
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The case of two sunlit tangent and two terminator solutions is an exception 
to the above test. There will be two segments sunlit separated by a dark segment. 
In most of these cases, Q 2 will fall in the dark segment. Then the total sunlit 
portion would still be the calculated terminator angle, 0 1 , minus the calculated 
tangent angle, 9 , that is (360 - d 1 ) - (360 -6) = 9 - 6 1 . 


VEHICLE 



Figure 4-Two Tangent and Two Terminator Plane Solutions 


E. Special Cases 


If two solutions are not found, the telescope cone is either completely sunlit 
or completely dark. In either case the complete cone is the same. Thus for some 
number k, \p + kS°| =R and the vector p + k S° will be sunlit if the cone is 
completely sunlit. Let p = (a, b, c), S° = (£ , m ,n ), then (a +k£ ) 2 + (b +km) 2 + 
(c +kn) 2 = R 2 which is a quadratic in k. Solve for the smallest k (in absolute 
value) and test if sunlit. If k is imaginary, S° does not intersect the central 
body and the cone is dark. (See Figure 5.) 


V. COMPUTER PROGRAM 

The method described has been programmed as a subroutine to the IMP Mis- 
sion Analysis Program,QUIMP, which is a modification of the Quick Look Mission 
Analysis Program, reference 1. The calling sequence is CALL TELANG (NOR) 
where NOR is the number of the central body (see QUIMP writeup). The following 
quantities are required from common, C(1000): 
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INTERSECTION OF 



C(58), C<59), C(60) 

C(21), .... C(28) 
C(721), .... C(744) 

C(138), . . . , C(140) 


The telescope-spin axis angles 

These must be input to QUIMP in OVRLAY: 

The radii of the ephemeris bodies. 

Position of the vehicle with respect to the ephemeris 
bodies. 

Spin axis vector. 


TELANG uses FNORM, DOT, and ADOT from the QUIMP program and 
subroutine QUART for solving the quartic equations. 


If two tangent solutions exist, TELANG computes the angle between them and 
prints this as the tangent solution. TELANG then tries to find two terminator 
solutions and prints their angle if it exists. If only one sunlit tangent solution 
exists, the tangent solution angle is printed as zero and the tangent solution is 
used with a terminator solution if it exists. 


VI. RESULTS 

Figures 6, 7, and 8 show some calculations done for the AIMP vehicle. Most 
cases were encountered, that is two tangent solutions, two terminator solutions, 
two of each, and one of each. The entire cone being sunlit was not encountered. 
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ANGLE THAT TELESCOPE SEES SUNLIT EARTH (degrees) 







180 



TIME FROM INJECTION (minutes) 





ANGLE THAT TELESCOPE SEES SUNLIT EARTH (degrees) 



11 




VII. CONCLUSIONS 


The described method has been successfully used as a subroutine in the 
QUIMP Program, and it could be used with any program which provides the 
necessary position and orientation information. 
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X. APPENDIX 


A. Calculations of Intersections (90 Degree Case) 


Let a x x + a 2 y + a 3 z + a 4 be the equation of the tangent or terminator plane 
as the case may be; bjX + b 2 y + b 3 z + b 4 =0 the equation of the telescope 
plane; and x 2 + y 2 + z 2 = R 2 the equation of the central body. Assume below 
that the coefficients divided by are not zero. If they are, solve for x or y instead 
of z . Then 


z = 


= (ajX + a 2 y + a 4 ) . 


( 1 ) 


Substituting (1) into the equation of the telescope plane, we get 


or 


bjX + b 2 y (a^ + a 2 y + a 4 ) + b 4 = 0, 


:jX + c 2 y + c 4 = 0 


( 2 ) 

(3) 


where the c's are obtained by collecting terms in (2). Solving for y , we have 


y (c, x + c ) 

C 2 


(4) 


and substituting this back into (1) , we obtain 


Z = - 


a 2 

a t x - — (CjX + c 4 ) + a 4 


(5) 


or 


z = djX + d 2 , 


( 6 ) 
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where the d’s are obtained by collecting terms in (5). Substituting for y and z 
from (4) and (6) into the equation of the central body, a quadratic equation in x 
is obtained. If there are any real solutions, use (4) and (6) to solve for y and z. 


B. Calculations of Intersections (non-90 degree case) 


Given the equation of the telescope cone (see p. 4) the equation of the tangent 
or terminator plane, i.e., 


b x x + b 2 y + b 3 z + b 4 = 0, (7) 

and the equation of the central body 

x 2 + y 2 + z 2 = R 2 , (8) 

we solve simultaneously. Again assume the coefficients divided by are not zero. 
Solving for z from (7), 

z = — g- (bjX + b 2 y + b 4 ), 
and substituting in (8), we obtain 


:,x 2 + c 2 y 2 + c 3 xy + c 4 x + c $ y + = 0 


(9) 


where 


b 2 b 2 2b. b 

= 1 + — , c = 1 + — , c, = 

b 2 2 b 2 3 b 2 


b? 5 b? 6 b 2 


3 3 

Consider (9) as a quadratic equation in y. Solving for y in terms of x; 


y = 


- (c 3 x + C 5 ) + /(C 3 X + C 5 ) 2 — 4c 2 (CjX 2 + c 4 x + c 6 ) 


2c n 


( 10 ) 


The - sign before the radical is superfluous in obtaining the final equation, but 
is necessary, once having solved for x , to obtain all intersections. Solving for 
z in terms of x , 


z = 




(-(c 3 x + c s )± J( 


C 3 X + C S ) " 4c 2 


(c x x : 


+ C 4 x + c 6 ) ) 



( 11 ) 
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Letting 


and 


y = djX + d 2 


+ d 3 / d 4 x 2 +d 5 x + d 6 


= e t x 


+ e o + e , 


/ 


e. x 

4 


e 5 X 


+ e * 


A 

1 





d S = e S = 2C 3 C 5- 4C 2 C 4 



- 4C 2 C 6 


e 


l 





b 

2b 


Expressions for y 2 and z 2 are needed for substituting in the telescope cone 
equation. 


where 


y 2 = fjX 2 + f 2 x + f 3 + (f 4 x + f 5 ) /d 4 x 2 + d $ x + d g 
z 2 = gjX 2 + g 2 x + g 3 + (g 4 x + g 5 ) /e 4 x 2 + e 5 x + e 6 


f l = d i + d 3 d 4 


f 4 = 2d i d 3 


+ e 3 e 4 


g 4 = 2e x e 3 


f 2 = 2d l d 2+ d 3 d 5 

f S= 2d 2 d 3 
g 2 = 2 ei e 2 + e 2 e s 

g 5 = 2e 2 e 3 


f 3 = d 2 + d 3 d 6 
g 3 = e 2 + e 3 e 6 


Now substituting in the equation for telescope cone 



where 


h x = a i + a 2 f l + a 3 *i + a 4 d l + a 5 e i + a 6 (d i e i + d 3 d 4 e 3 ) 

h 2 = a 2 f 2 + a 3 g 2 + a 4 d 2 + a 5 e 2 + 3 6 (d l C 2 + d 2 e i + d 3 e 3 d 5> + 3 7 + Vl + 3 9 e ! 

h 3 = a 2 f 3 + a 3«3 + a 6 (d 2 e 2 + d 3 e 3 d 6> + 3 8 d 2 + 3 9 e 2 + 3 10 
h 4 = a 2 f 4 + 3 3 ®4 + a 4 d 3 + 3 S e 3 + 3 6 (d l e 3 + d 3 e i> 

h 5 = a 2 f 5 + a 3 g S + 3 6 (d 2 e 3 + d 3 e 2 ) + 3 8 d 3 + a 9 e 3 

Squaring both sides above results in a quartic equation in x 

PjX 4 + p 2 x 3 + p 3 x 2 + p 4 x + p s = 0 
where. 

Pi = h 2 - h 2 d 4 P 2 = 2h 2 h! - (h 2 d g + 2h 4 h 5 d 4 ) 

p 3 = 2h!h 3 + h 2 - (h 2 d 6 + 2h 4 h 5 d s + h 2 d 4 ) p 4 = 2h 2 h 3 - (2h 4 h s d 6 + h 2 d s ) 

Ps = h 3- h s d 6 

The roots of this quartic, and equations (10) and (11) give the solutions. 
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